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ABSTRACT 

A model of supernova-driven chemical evolution of the Galactic halo, recently pro- 
posed by Tsujimoto, Shigeyama, & Yoshii (1999, ApJL, 519, 63), is extended in order to 
investigate the evolution of light elements such as Be and B (BeB), which are produced 
mainly through spallative reactions with Galactic cosmic rays (GCRs). In this model 
each supernova (SN) sweeps up the surrounding interstellar gas into a dense shell and 
directly enriches it with ejecta which consist of heavy elements produced in each Type 
II supernova with different progenitor masses. We propose a two-component source for 
GCRs such that both interstellar gas and fresh SN ejecta engulfed in the shell are acceler- 
ated by the shock wave. The released GCRs travel much faster than the expansion of the 
shell, and thus produce the BeB elements far outside the shell which will be incorporated 
in subsequent formation of shells arising from later SNe. As a result, stars formed from 
coeval shells are predicted to show a large scatter in the abundance of heavy elements, 
while exhibiting BeB abundances similar to that in the gas, with no appreciable scatter. 
This indicates that, contrary to heavy elements, stellar BeB abundances might be used 
as a good age indicator in the inhomogeneous Galactic halo. 

The production of BeB at early epochs is dominated by the primary process through 
spallation of heavy GCRs, though being a minor component in the bulk of the GCR 
composition at present. We have calculated the frequency distribution of long-lived stars 
in the log(BeB/H)-[Fe/H] plane and find that the contour of constant frequency covering 
a range of —3 <[Fe/H]< —1 in this plane is consistent with the observed linear trend 
between BeB and Fe. We show from our calculations that there is an intriguing possibility 
of distinguishing between standard and non-standard Big-Bang nucleosynthesis models 
if BeB abundances in several hundred halo stars are observed in the future. 
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1. INTRODUCTION 

The study of elemental abundance patterns in very 
metal-poor halo stars has suggested that stars with 
—4 <[Fe/H]< —2.5 were born from individual su- 
pernova remnants (SNR), reflecting a different nu- 
cleosynthesis yield for each Type II supernova (SN 
II) with different progenitor mass (Shigeyama & Tsu- 
jimoto 1998, hereafter ST98; Ishimaru & Wanajo 
1999). This indicates that previous models of chemi- 
cal evolution in which stars are assumed to form from 
well-mixed, uniform gas cannot apply to the Galactic 
halo; as a result a unique age versus metallicity re- 
lation is not expected to hold for stars with [Fe/H]< 
—2.5. Tsujimoto, Shigeyama, & Yoshii (1999, here- 
after TSY) have recently presented a new model for 
chemical evolution in the early Galactic halo, explic- 
itly taking into account the expected inhomogeneous 
nature of the gas, and showed that their model is 
useful in interpreting the stellar abundances of heavy 
elements originating from SNe II. 

The light elements such as Be and B (BeB) in halo 
stars, which are the focus of this paper, have usu- 
ally been considered to originate from energetic pro- 
tons and a-particles in Galactic cosmic rays (GCRs) 
colliding with heavy nuclei of C, N, and O (CNO) 
in the interstellar medium (ISM) (e.g., Walker et al. 
1985). While this secondary process naturally pre- 
dicts a quadratic dependence of BeB abundance on 
metallicity, newly obtained data of stellar BeB abun- 
dances have, to the contrary, established a linear de- 
pendence (Boesgaard et al. 1999 for Be; Duncan et 
al. 1997 and Primas et al. 1999 for B). 

To account for this discrepancy, many authors have 
considered possible primary processes for producing 
the BeB elements (e.g., Olive et al. 1994; Casse et 
al. 1995; Vangioni-Flam et al. 1996; Lingenfclter et 
al. 1998; Vangioni-Flam et al. 1998; Higdon et al. 
1998). In particular, Duncan et al. (1992) suggested 
that GCR CNO collide with protons and a-particles 
in the ISM to produce BeB, which leads to their linear 
dependence on metallicity (see also Yoshii, Kajino, & 
Ryan 1997; Ramaty et al. 1997). From a reanalysis of 
the elemental abundances of present-day GCRs and 
an examination of the physics of the SNR shock wave 
acceleration process, Meyer & Ellison (1999, hereafter 
ME99) and Ellison & Meyer (1999) have recently ar- 
gued that most GCRs were accelerated out of the 
ISM or circumstellar material instead of SN ejecta, 
although they suggested that, owing to the lack of 



ISM CNO at early epochs, GCR CNO out of fresh 
SN ejecta might play a role in producing the BeB 
elements. 

In this paper, by incorporating both primary and 
secondary processes of GCR spallation of heavy ele- 
ments, we can, for the first time, construct a model 
for the evolution of light elements which is fully con- 
sistent with the inhomogeneous chemical evolution of 
the early Galactic halo. The novelty of the result 
based on this model is that the frequency distribution 
of metal-poor stars in the log(BeB/H)-[Fe/H] plane 
is predictable, along with the expected intrinsic abun- 
dance scatter. 

2. FORMULATION 

2.1. Inhomogeneous Chemical Evolution 

We assume that the star-forming process is con- 
fined in separate clouds, which make up the entire 
halo of the early galaxy, and that the GCRs arise 
from multiple SN explosions therein. In our model a 
certain fraction ini of cloud mass has turned into Pop 
III stars at time t = and the massive stars among 
them explode as Pop III SNe to initiate the chemical 
evolution in a given cloud. All stars of subsequent 
generations are assumed to form from the SNR shells 
behind the radiative shock front, which is supported 
by recent observations of the abundance patterns in 
metal-poor stars (ST98). The mass fraction of each 
shell turning into stars is taken as a constant, denoted 
by e. Following TSY, provided xm < 1 x 1CT 2 , a value 
of e = 4.3 x 1CP 3 is used in the star formation rate 
given by 

M*(i)= ( dmeM sh (m,t)^^M*(t-T{m)), (1) 

Jmax(m t ,m SNi i) m 

where t(to) denotes the lifetime of a star with mass 
m, and rrij is the stellar mass for which 7"(m) = t. The 
initial mass function, 4>(m), having a Salpeter slope 
index of —1.35, is adopted with upper and lower mass 
limits of m u — 50Mq and to; = 0.05Mq. A lower 
mass limit of stars that explode as SNe is taken to 
be TOsn,2 = IOMq. The mass of the shell is given by 
M sh (TO,i) = M ej (ra) + M SVf (m,t), where M ej (m) is 
the mass of the SN ejecta and M sw (m, t) is the mass 
of the gas swept up by the SNR. Since M SVJ (m,t) is 
insensitive to m and t, we here adopt M sw (m,t) = 
6.5xIO 4 M (ST98; TSY). 

The mass of the gas, M g (t), and the abundance 
of the j-th heavy element in the gas, Zj ig (t), change 



2 



with time according to eqs.(2) and (3) in TSY, re- 
spectively. Because of the inhomogeneous mixing of 
chemical compositions in the dense shells and the 
ISM, the abundance of the j-th heavy element in 
stars born at time t in an SNR shell depends on the 
mass m of SN progenitor, and is essentially given by 
Z^(m,t) = [M Zj (m) + Z jtS (t)M BVf (m,t)]/M Bh (m,t) 
(see TSY). 

Since the light elements (LiBeB) are easily de- 
stroyed by stellar nuclear reactions at temperatures 
above a few 10 6 K, the abundance of the L element in 
the gas Zi, lg (t) changes through destruction by stellar 
lock-up and production by GCR reactions according 
to: 



d(Z Lyg M g 
dt 



4>(m) 



— / dmZL,*(m,t)eM s h(m,t) 

'max(rot,msN,i) m 



xM»(t - r(m)) 



]T { F% ) Zj,g (t) [A L /Aj ) 

=pa,j=CNO 



{a^X^AL/A^M^t), 



where A, is the mass number of the i-th element and 
Xi is the abundance of hydrogen or helium. We define 
(4^) ee J^a^F^E^S^dE, where a^E) 
is the cross section for the process of GCR projec- 
tile i impinging on the ISM target j to produce the 
L element. Sl(E) gives the retention fraction of L 
products that can survive to be thermalizcd into the 
ISM, and F l (E,t) is the time-dependent flux of GCR 
projectile i, which is modeled in §2.2. The stellar L 
abundance is equal to that in the shell and is given 

by 

Zi,*(m, t) = [Z L ,g(t)AW(m, t) + 

i=pa,j=C~NO 

{ { a ij ^ ) Zj,g (t) (Al /Aj ) + (af^XiitXAL/Ai)} 

x M sh (m,t)AT]/M sh (m,t) , (3) 

where AT = 10 6 yr is a typical diffusion time of a 
SNR shell (Shigeyama 1999, private communication). 
The first term in the numerator represents the mass of 
L element swept up in the shell and the second term 
represents that mass produced during AT by GCRs 
originating from all the SNe that explode at time t. 

2.2. A Model of Time-Dependent GCR Flux 

The energy spectrum of the flux of propagated 
GCR projectiles, denoted by Fi(E,t) in units of 
particles s _1 cm _2 (MeV per nuclcon) -1 , is expresssed 
as Fi(E,t) = M£)]"V£ *(£',*) expHi^OE') - 



Ri(E))/A]dE' (Meneguzzi et al. 1971; Prantzos et 
al. 1993), where qi(E,t) is the source spectrum in 
particles s -1 g -1 (MeV per nucleon) -1 , A is the effec- 
tive escape length in gem -2 , uJi(E) is the ionization 
energy losses in MeV(g cm -2 ) -1 through a hydrogen- 
helium plasma with X=0.75 and Y=0.25, and Ri(E) 
is the ionization loss length in g cm -2 (Northcliffe & 
Schilling 1970). The effective escape length is given 
by A -1 = A -1 + A -1 where A c is the loss length of en- 
ergetic particles against escape from the Galactic halo 
and A n is that against nuclear destruction and decay. 
A e is proportional to the scale height if in a cylin- 
drical collapse of the halo (Prantzos et al. 1993) and 
approximately changes with time as A c oc 1 — (at) 2 , 
where a -1 is the timescale of the collapse and a value 
of 1 Gyr is used below. Although A e would be huge 
in the past, the maximum value of A is bounded by 
A n ~ 200 g cm -2 for protons and 20 g cm -2 for a 
particles (Malaney & Butler 1993). 

We propose a new scenario in which GCRs origi- 
nate from a combination of injected SN ejecta (pri- 
mary origin) and swept-up ISM (secondary origin) in 
the SNR shell, which are then accelerated through 
the scattering back and forth across the shock. Their 
source spectrum qt(E, t) is assumed to have the same 
form, with its intensity proportional to the ejection 
rate of energetic particles associated with SNe, that 
is: 



qi {E,t)<x{E + E ) 



-2.6 



dm 

max(m t , msN.l) 



x{M Zi (m)+Z itg (t)f cr M sw (rn,t)}^-tM*{t - r(m)), 

Aim 

(4) 

where Eq = 931MeV, M^fm) is the mass of the i- 
th heavy element synthesized and ejected from a star 
with mass m, and f cr is the fractional mass of the gas 
in the shell swept up while the SN explosion is able 
to accelerate ISM particles. If the acceleration were 
effective during the Sedov phase of the adiabatic shell 
expansion, a typical value of f cr would be 0.01 — 0.1 
(Taylor 1950; Sedov 1959). The absolute value of q t 
is chosen so that the predicted Be abundance from Fi 
should agree with the observed abundance at [Fe/H] = 
—3. The composition of a particles relative to protons 
in GCRs is assumed to remain equal to a solar value of 
F a /F p = 0.135 (Lund 1989), but the ratio F CNO /F p 
is calculated by using eq.(U). 
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3. RESULTS AND DISCUSSION 



We have calculated the evolution of light and heavy 
elements in an inhomogeneous Galactic halo using 
the formulae of the preceding section. The val- 
ues of the free parameters chosen in our model are: 
x m = 1CT 6 , e = 4.3 x 1CT 3 , and f cr = 0.01. The ini- 
tial abundances of heavy elements are set to zero, but 
those of BcB elements are taken as equal to the pre- 
dicted primordial abundances of log(Bc/H)=— 17.9 
and log(B/H)=— 16.9, based on the standard Big 
Bang nucleosynthesis (SBBN) model by Thomas et 
al. (1994). 

Our results are shown in the color-coded frequency 
distribution of long-lived stars with m < 1-Mq in the 
log(Be/H)-[Fe/H] plane (Fig.l) and in the log(B/H)- 
[Fc/H] plane (Fig. 2). Both sets of frequency dis- 
tributions have been convolved with a typical error 
a = 0.15 dex in the abundance measurements in order 
to directly compare with the Keck data for Be (Boes- 
gaard et al. 1999) and the HST data for B (Duncan 
et al. 1997; Primas et al. 1999). 

We define a probability density of finding one halo 
star within a unit area of A[Fe/H]=0.1 x Alog(BeB/H) 
=0.1, normalized to unity when integrated over the 
entire area. The five contour lines shown in Figs.l 
and 2 are, from the inside out, of constant proba- 
bility density 10~ 2 , 10" 3 , 10~ 4 , 10~ 6 , and 10~ 8 ; the 
total probability integrated over the area inside each 
contour is 0.57, 0.95, 0.991, 0.9999, and 0.9999997, re- 
spectively. Inspection of these figures shows that the 
area of 0.95 probability is in reasonable agreement 
with the distribution of observed data for both Be 
and B, and that, except for the presumed existence of 
metal-free Pop III stars of the first generation, there 
should exist no stars below [Fe/H]« —5 because the 
SNR shells at each star formation site are enriched 
by the heavy-element ejecta from Pop III SNe. We 
note, however, that this value on the lower bound in 
the stellar metallicity is subject to the nucleosynthesis 
calculations adopted for SNe II. 

The pattern of predicted abundance scatter in 
the log(BeB/H)-[Fe/H] plane changes distinctly from 
lower to higher metallicity regions. This transition oc- 
curs at [Fc/H]~ —2.5, below which the heavy-element 
mass of SN ejecta exceeds the heavy-element mass in 
the gas swept up in the SNR shell (ST98). Therefore, 
in the lower metallicity region ([Fe/H]< —2.5), stars 
formed at the same epoch from different SNR shells 
will exhibit a range of metallicities, reflecting a differ- 
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Fig. 1. — Color-coded frequency distribution of long- 
lived stars in the log(Bc/H)-[Fe/H] plane, convolved 
with a Gaussian having a = 0.15 dex for both 
log(Be/H) and [Fe/H]. The primorial Be abundance 
is taken from the SBBN model, and is indicated by 
the dashed line. The five contour lines, from the in- 
side to the outside, correspond to those of constant 
probability density 10" 2 , 10" 3 , 10" 4 , 10" 6 , and 10" 8 
in unit area of A [Fe/H] =0.1 x Alog(Be/H)=0.1. The 
thick line shows the log(Be/H)-[Fe/H] relation for the 
gas for our standard choice of f cr = 0.01; the lines in 
the inset show those for various values of f cr . The 
crosses represent the data with observational errors 
taken from Boesgaard et al. (1999). 

ent nucleosynthesis yield for each SN II with differ- 
ent progenitor mass. On the other hand, the GCRs, 
which are accelerated either directly from each SN II 
or from swept-up ISM, can escape from the SNR shell 
at high speed without producing BeB through spalla- 
tive reactions in the shell. As a result, most of BcB 
are produced in the ambient ISM far outside the shell, 
although the GCRs coming from outside impinge on 
targets in the shell to produce a small amount of BcB 
there. This accounts for the prediction that stars 
with [Fe/H]< —2.5 show a large scatter in the heavy- 
element abundance (along the horizontal axis), while 
exhibiting only a small scatter in the BeB abundance 
(along the vertical axis). It is therefore important to 
note that purely GCR-induced products like BeB can 
still be used as an good age indicator at early epochs, 
although the so-called age versus metallicity relation 
is broken. 

The predicted scatter in [Fe/H] becomes smaller 
towards higher [Fe/H], so that the frequency distribu- 
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Fig. 2. — Color-coded frequency distribution of long- 
lived stars in the log(B/H)-[Fe/H] plane, as in Fig.l 
but for B ( 10 B + n B). Crosses represent the data 
taken from Duncan et al. (1997) and Primas et al. 
(1999) 

tion of stars in the log(BeB)-[Fe/H] plane converges to 
a linear relation BeB/HocFe/H, as is observed over a 
range of —3 <[Fe/H]< —1. This linear trend emerges 
when the primary GCR processes produce BeB. No- 
tably, as suggested by ME99, the primary processes 
by GCR CNO impinging on ISM pa, which are over- 
whelmed by the secondary processes at later epochs 
(including the present), become dominant at early 
epochs because of the lack of ISM CNO. This does 
not mean that most GCRs are of the primary origin. 
The mass ratio of primary to secondary GCRs in our 
model is M c j{m)/(f cr M sw (m,t)) w 10/(0.01 x 6.5 x 
10 4 ) ~ 0.02 (cf. eq.([|)), which is not inconsistent with 
the recent results by ME99 that the GCRs originate 
mostly from the ISM. This ratio could also be argued 
by energetics between the kinetic energy released from 
SN and the input energy onto GCRs(Ramaty et al. 
1997), which will be discussed in detail in our forth- 
coming paper. 

As [Fe/H] increases, the secondary GCR processes 
become dominant and a quadratic relation BcB/Hoc 
(Fe/H) 2 results. This changeover from a predicted 
linear to quadratic trend occurs at lower metallicity 
for larger f cr , as shown in the inset of Fig.l. If a 
larger f cr is used, in a typical range of 0.01 — 0.1 
(Taylor 1950; Sedov 1959), the predicted relation de- 
viates upwards from the observed linear trend. This 
deviation is partly offset by the effect of gas outflow 
from the cloud, which is known to flatten the slope of 



the log(BeB/H)-[Fe/H] relation in the late phase of 
chemical evolution (e.g., Yoshii et al. 1997). In this 
paper we have taken into account this outflow only 
approximately by using a lower value of f cr = 0.01 
as our standard. More BeB data for stars in the in- 
termediate metallicity region —2 <[Fe/H]< —1 are 
highly desirable to constrain f cr , that is, the acceler- 
ation mechanism of GCRs as well as the rate of gas 
outflow from the cloud. 

So far we have used the standard BBN model to 
set our assumed initial abundances of BeB elements, 
although their abundances are strongly dependent on 
the BBN model which is chosen. For example, if the 
density distribution in the early universe is not uni- 
form as in the inhomogeneous BBN (IBBN) model, 
more LiBeB arc synthesized in neutron-rich regions 
(Kajino & Boyd 1990). In order to see how a non- 
standard BBN model changes the result, we calcu- 
late the frequency distribution of long-lived stars in 
the log(Be/H)-[Fe/H] plane using an optimal initial 
abundance of log(Be/H)=— 14.3 taken from the IBBN 
model (Orito et al. 1997). The resulting distribu- 
tion for the IBBN model is compared with the SBBN 
model in Fig. 3, where the contour lines of constant 
probability density 10~ 4 and 10~ 5 in unit area of 
A[Fe/H]=0.1xAlog(Be/H)=0.1 are shown. Since the 
total probability integrated over the area between the 
inner and outer contours is less than 0.01, we need 
to measure BeB abundances in several hundred halo 
stars in order to distinguish the IBBN from the SBBN 
on a statistical basis. A more direct way to distiguish 
these alternative hypotheses is to identify true Pop III 
stars, the atmospheres of which consists of primeval 
gas. According to TSY, at least one halo star out of 
10 3-4 could be such a metal-free Pop III star. Such a 
large sample already exists (e.g., Beers 1999) and the 
discovery of Pop III stars may not be an unrealistic 
wish in the near future. 

4. SUMMARY 

We have presented a general approach for studying 
the GCR-induced evolution of the L elements (LiBeB) 
in a manner consistent with an inhomogeneous chem- 
ical evolution of the Galactic halo. By adoption of 
a two-component source of GCRs, our model results 
include: (1) a prediction of the intrinsic scatter in 
BeB and [Fe/H] abundances within the model, (2) a 
successful prediction of the observed linear trend be- 
tween BeB and Fe, (3) a proposal for using BeB as 
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Fig. 3. — Frequency distribution of long-lived stars 
in the log(Be/H)-[Fe/H] plane for two different pri- 
mordial Be abundances predicted in the SBBN and 
1BBN models. The contour lines, from the in- 
side to the outside, correspond to those of constant 
probability density 10~ 4 and 10~ 5 in unit area of 
A[Fe/H]=0.1xAlog(Be/H)=0.1 for the SBBN (outer 
boundaries of two regions painted by different colors) 
and for the IBBN (dashed lines). The log(Be/H)- 
[Fc/H] relation for the gas is shown by the thick line 
for the SBBN and by the dashed line for the IBBN. 
The crosses represent the data with observational er- 
rors taken from Boesgaard et al. (1999). 

a cosmic clock, as an alternative to [Fe/H], and (4) 
a method for possibly constraining the BBN model 
from future observations of metal-poor stars. 
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